To APPEARIN The Astrophysical Journal, Letters. 
Preprint typeset using LATgX style emulateapj v. 21/08/00 



O 
O 

o 

O 

Q 

o 



FIRST RESULTS FROM THE ANGLO- AUSTRALIAN PLANET SEARCH 

CANDIDATE AND A 51 PEG-LIKE PLANET 1 



A BROWN DWARF 



C.G. 



Tinney 2 , R. Paul Butler 3,2 , Geoffrey W. Marcy 4,5 , Hugh R.A. Jones 6 , Alan J. 
Penny 7 , Steven S. Vogt 8 , Kevin Apps 9 , Gregory. W. Henry 1 " 

To appearin The Astrophysical Journal, Letters. 

ABSTRACT 

We report results from the Anglo- Australian Planet Search - a survey for planets around 200 solar-type 
stars in the southern hemisphere, which is being carried out on the 3.9m Anglo-Australian Telescope. 
Limiting Doppler precisions of 3ms _1 have been demonstrated from the first 2.5 years of operation, 
making this the highest precision planet search in the southern hemisphere. From these data we report 
results for two new sub-stellar detections. The first is a "51Peg"-like planet around the star HD 179949 
with Msini = 0.84 Mjtjp- Photometric study reveals this is not a transiting system. The second is a 
brown dwarf or very low-mass star companion to HD 164427 in an eccentric orbit with M sin i = 46 Mjup ■ 
Hipparcos data indicate this latter object is unlikely to have a mass greater than O.18M0. 
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pixel CCD. This lumogen-coated CCD (denoted MITLL2) 
underwent a serious failure at the end of 1999, which re- 
quired its read-out to be switched to the second of its 
working amplifiers in 2000 January (following which it was 
denoted MITLL2a). No change in the operation of the de- 
tector (apart from a flipped read-out format) has been 
detected by our programme. 

Our target sample of «200 stars with 6 < —20°, in- 
cludes F,G and K V-IV stars with V<7.5 and MV stars 
with V<11.5. Where age/activity information is available 
from S or Rhk indices (see e.g., Henry et al. (1996)), 
target stars are required to have ages greater than 3 Gyr. 
Our first observing run was in 1998 January, and the last 
run for which observations are reported here was in 2000 
November. The observing and data processing procedure 
follows that described in Butler et al. (1996), and is 
described in detail in Paper II. In particular, Paper II 
presents velocities for a number of stable stars illustrating 
that we reach a velocity precision floor of 3ms _1 for the 
bright stars included in our programme. One of the stars 
presented here (HD 179949) has also had confirmatory ob- 
servations acquired at three epochs in 2000 September 5-8, 
as part of the Keck planet search (Vogt et al. 2000). 

3. CHARACTERISTICS OF HD 179949 & HD 164427 

HD 179949 (HR 7291, HIP 94645, GJ749) is an F8V star, 
which shows moderate rotation with a measured v sin i = 
6.3±0.9kms _1 (Groot et al. 1996). Its Hipparcos par- 
allax puts it at 27.0±0.5pc, with M y =4.09±0.04 (ESA 

1 Based on observations obtained at the Anglo— Australian Telescope, Siding Spring, Australia; and the W.M. Keck Observatory, which is 
operated jointly by the University of California and the California Institute of Technology. 
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1. INTRODUCTION 

Since the discovery of the first extra-solar planet by 
Doppler velocity techniques in 1994 by Mayor & Queloz 
(1995), planetary detections have been dominated by 
northern hemisphere search programmes - most prolifi- 
cally by the precision velocity programmes at Lick (e.g. 
Butler et al. (1996)) and Keck (e.g. Vogt et al. (2000)), 
but also complemented by lower precision programmes 
at OHP (Baranne et al. 1996), McDonald Observatory 
(Cochran et al. 1997), AFOE (Noyes et al. 1997), and 
programmes at La Silla (Kiirster et al. 2000; Queloz et al. 
2000). Of these programmes, only the latter have access 
to the sky south of ~ —20°, and these achieve precisions 
of ~10ms _1 . In 1998, therefore, the Anglo-Australian 
Planet Search was begun to complete the all-sky cover- 
age of the brightest stars at precisions reaching 3ms" 1 . 
In this paper we present some first results from this pro- 
gramme. A companion paper (Butler et al. 2000, here- 
after Paper II) presents results for a further two new plan- 
ets, along with a detailed description of our observational 
programme. 

2. THE ANGLO-AUSTRALIAN PLANET SEARCH 

The Anglo- Australian Planet Search is being carried out 
on the 3.92m Anglo- Australian Telescope (AAT), using 
the University College of London Echelle Spectrograph 
(UCLES) and an I2 absorption cell. UCLES is operated in 
its Sllinesmm" 1 mode with an MIT/LL 2048x4096 15^m 
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1997). There is no published evidence to indicate that 
HD 179949 is a binary. The photometry of Eggen (1998) 
indicates HD 179949 has roughly solar mctallicity with 
[Fe/H]=+0.02±0.1. The ubvy calibration of Vogt et al. 
(2000) and published photometry would indicate a more 
metal-rich [Fe/H]=+0.22±0.07. The latter is supported by 
Twarog & Anthony- Twarog (1995, Fig. 8) which indicates 
HD 179949 is as metal-rich, if not more so, than the major- 
ity of disk main-sequence stars. HD 179949 has a moder- 
ate X-ray luminosity as determined by ROSAT's PSPC - 
Hiinsch et al. (1998) determine a Lx=41.0xl0 27 ergs -1 , 
while in an independent analysis Piters et al. (1998) deter- 
mine Lx/Lfc o ;=(9.7±3.3)xl0~ 6 . Both are some 10 times 
higher than the equivalent quantity in the quiet Sun. The 
spectra acquired at Keck have been used to derive a R' hk 
index of -4.72 for this star, using the same procedure as 
Vogt et al. (2000). Figure 1 compares the Ca II H line 
with that seen in the Sun. The mass of HD 179949 is esti- 
mated to be 1.24±0.1M Q , based on interpolation between 
evolutionary tracks of Fuhrman et al. (1997, 1998). 

HD 164427 (HIP 88531, G1700A) is an inactive GOV 
(Houck & Cowley 1975) star, with a R' hk index of - 
4.95 (Henry et al. 1996). Its Hipparcos parallax puts 
it at a distance of 39.1±1.4pc, and M y =3.91±0.8 (ESA 
1997) making it somewhat over-luminous for its spec- 
tral type. Indeed Evans et al (1964) classified it as a 
sub-giant, with luminosity class IV. The uvby calibra- 
tion of Vogt et al. (2000) suggests a metallicity of 
[Fe/H]=+0.11±0.07. HD 164427 was catalogued as a very 
wide binary (G1700AB) by Gliese (1969) (though both 
have since dropped out of the nearby star sample (Gliese & 
Jahreiss 1991) due to better parallaxes placing them out- 
side the 25 pc limit). The binary identification is actually 
due to Luyten (1957) who catalogued this pair as LTT 7172 
and LTT 7173 (respectively) with common proper motion 
and a 28" separation at a position angle of 336°. There 
is no subsequent astrometry of this pair in the literature. 
Examination of Digital Sky Survey scans of a UKST plate 
from 1992.6 shows no evidence for an object with the mag- 
nitude difference indicated by Luyten (Am = 7.2) at this 
separation relative to HD 164427, leading us to conclude 
it may not be a common proper motion pair. Even if 
the system were a binary, it would be a very wide system 
(1090 au) - so wide as to be irrelevant for the purposes of 
high precision doppler velocities. The mass of HD 164427 
is estimated to be 1.05±0.1M©. 

Both stars were seen to be photometrically stable over 
the life of the Hipparcos mission at a 95% confidence level 
of <0.015 magnitudes (ESA 1997). 

4. RADIAL VELOCITY OBSERVATIONS AND ORBITAL 
SOLUTIONS 



effects of photon-counting uncertainties, residual errors in 
the spectrograph PSF model and variation in the underly- 
ing spectrum between the template and "iodine" epochs. 
Only observations where the uncertainty is less than twice 
the median uncertainty are listed. These data are shown 
in Figure 2 phased with a period of 3.093 days. The fig- 
ure shows the best-fit Kcplcrian model for the data, with 
the resultant orbital parameters listed in Table 2. The 
best-fit period from Keplerian fitting is the same as the 
period found from the peak of the Scargle periodogram, 
with a false alarm probability less than 0.001. Both A AT 
(dots) and Keck (squares) data are shown, with the Keck 
data completely confirming the orbital fits derived from 
the AAT data alone. Like the previously observed giant 
planets in 3 to 5 day orbits, HD 179949 appears to be in 
a circular orbit. Table 2 includes the best-fit orbital so- 
lutions both for the case of a circular orbit and floating 
eccentricity. The resulting minimum companion mass is 
0.84Mjup, with an orbital semi-major axis of 0.045 au. 

The residuals about the fit are at the level of "jitter" ex- 
pected in an F8V star with HD 179949's level of rotation 
(e.Skms- 1 ) and activity (i^ fe =-4.72). "Jitter" here is 
used to refer to the scatter in the observed velocity about 
a mean vauc in systems observed over the long-term to 
have no Keplerian Doppler shifts. It is thought to be the 
combined effect of surface inhomogeneties, stellar activ- 
ity and stellar rotation. Saar et al. (1998) have studied 
the correlations between "jitter" , and stellar rotation and 
activity in the Lick precision velocity programme. They 
indicate we would expect HD 179949 to show «10ms -1 
"jitter" due to activity, spots, etc. We measure a scatter 
of 10 ms -1 about our Keplerian fit, which is consistent 
with this expectation. 

The twenty-seven observations of HD 166427 are listed 
in Table 3, and they are shown in Figure 3 along with a Ke- 
plerian fit to the data with the orbital parameters listed in 
Table 2. The rms scatter about this fit is 7.7 m s -1 , slightly 
larger than the median uncertainty from the fitting process 
of 5.6 ms -1 . This is larger than the systematic precision 
limit of the Anglo- Australian Planet Search, as the ear- 
liest observations of this object were performed in very 
poor conditions. As soon as it was realised HD 164427 
was a velocity variable with large amplitude, it was "pro- 
rated" within the programme to only short exposures (and 
so lower precision than the rest of the programme), though 
still more than sufficient to precisely determine the com- 
panion's orbit. The resultant minimum companion mass 
is 46Mjup, or 0.043M Q , and the orbital semi-major axis 
is 0.46 au. 

5. DISCUSSION 

5.1. HD 179949 



Twenty-three observations of HD 179949 are listed in Ten other "51 Peg" -like planets are currently known 

Table 1, where the column labelled "Uncert." is the veloc- with periods under 10 days, and eccentricities less than 

ity uncertainty produced by the least-squares fitting pro- 0.1 11 . The minimum masses of these planets span the 

cess, which simultaneously determines the Doppler shift range 0.24-3.87 Mjup, with orbital semi-major axes in 

and the spectrograph point-spread function (PSF), given the range 0.038-0.059 au. In all respects, HD 179949 falls 

an iodine absorption spectrum, an "iodine free" template squarely into this class of objects. Even its above-solar 

spectrum of the object, and an "iodine" spectrum of the metallicity is similar to that seen in the other "51Peg"- 

object (Butler et al. 1996). This uncertainty includes the like planets (e.g., Queloz et al. (2000), Gonzalez (2000)). 

11 In the absence of a formal classification scheme, we consider HD 83443, HD 168746, HD 46375, HD 75289, BD-10 3166, HD 187123, HD 209458, 
u And b, tBoo to be "51 Pcg"-like, as well as 51 Peg itself. 
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5.2. Transit Search in HD 179949 

As with all the other "51 Peg" -like planets, HD 179949 
is an excellent candidate for a transit search. The Hip- 
parcos magnitude, color index, and parallax of the star 
imply a radius of about 1.24 R Q , which, combined with 
the orbital radius of 0.045 au, places the probability of 
observable transits at 12.8%. Successful transit observa- 
tions (already carried out for HD 209458 by Charbonneau 
et al. (2000) and Henry et al. (2000b)) not only tightly 
constrain orbital inclination, but provide planetary radii 
as well. 

To search for transits in HD 179949, we obtained photo- 
metric observations between 2000 September 21 and Oc- 
tober 1 with the T8 0.80 m Automatic Photoelectric Tele- 
scope (APT) at Fairborn Observatory in southern Ari- 
zona 12 . The instrumentation and techniques used for this 
transit search are described in Henry (1999) and Henry 
et al. (2000a) . The observations were made differentially 
with respect to the comparison star HD 178075 (V = 6.30, 
B-V = 0.01, B9.5V). The precision of an individual 
observation is approximately 0.004 mag, somewhat worse 
than the typical precision of 0.001 mag obtainable with the 
APT due to observations being at air masses between two 
and three. 

No transit events were detected at the times predicted 
by the orbital parameters of Table 2. For an assumed plan- 
etary radius of 1.4 Rjup (the measured radius of the plane- 
tary companion of HD 209458 (Henry et al. 2000b; Jha et 
al. 2000)), the predicted transit depth is 0.014mag. The 
mean of our 16 photometric observations of HD 179949 
taken within the predicted transit window is 0.1058 ± 
0.0010 mag; the mean of the 38 observations outside the 
transit window is 0.1044 ± 0.0007 mag. Thus, the two 
means agree within their respective errors and transits 
deeper than about 0.001 mag are ruled out. This non- 
detection of transits limits the orbital inclination i to less 
than 83°, and sini to less than 0.992. Furthermore, the 
photometric observations place a limit of about 0.001 mag 
on any photometric variability of HD 179949 on the radial 
velocity period. This confirms that neither starspots nor 
stellar pulsations can be the cause of the radial velocity 
variations (see Henry et al. (2000a)) and so strongly sup- 
ports the existence of the planet even in the absense of 
transits. 

5.3. HD 16U27 

Recent work on the mass function of brown dwarfs in 
star clusters (Luhman et al. 2000; Lucas & Roche 2000; 
Zapatero Osorio et al. 2000) has failed to show that 
star formation processes cannot form objects below the 
13Mjup deuterium burning limit. (Though the same re- 
sults can be interpreted - less straightforwardly - as im- 
plying the existence of free-floating sub-13Mjup objects 
formed by dynamical evolution ejecting objects formed by 
planetary processes.) Certainly they provide no evidence 
that the deuterium burning limit is a useful demarcation 
boundary between star- and planet-formation processes. 
Difficult as it may be, only an understanding of forma- 
tion mechanisms (and possible subsequent dynamical evo- 
lution) can allow us to determine whether objects in the 
5-20 Mjup range are planets or brown dwarfs. Nonethe- 

12 These observations are available at http : //schwab. tsuniv . edu/t; 



less, the minimum mass of the companion to HD 164427 
clearly lies in a range unlikely to be due to planetary for- 
mation processes - it is a brown dwarf or a star, not a 
planet. 

Hipparcos astromctry, however, allows us to place some 
limits on the orbital inclination and mass for HD 164427's 
companion. The astrometric solution for HD 164427 shows 
a good fit for the parallax and proper motion of a single 
star, with 1-a residuals of «3mas. The uncertainty on the 
derived parallax is 0.9 mas. Based on this we can assume 
an upper limit to the astrometric perturbation (ao) of 
HD 164427 due to an unseen companion of ~ 2 mas. Then 
from the relationship between the spectroscopic and as- 
trometric orbital elements (e.g. Halbwachs et al. (2000), 
equation 1) we conclude that sinz > 0.24, or equivalcntly 
that the companion mass is less than 190Mjup or 0.18 M Q , 
making it either a brown dwarf, or a very low-mass star. 
With an orbital semi-major axis of 0.46 au at 39 pc, such a 
companion will have an maximum apparent separation of 
11.8 mas. Table 4 compares the absolute K magnitudes of 
HD 164427, with predicted magnitudes of the companion 
in the range of allowed masses. Only at the very mas- 
sive end of this range is the companion likely to be de- 
tectable by either ground-based adaptive optics or space- 
based imaging. 

Nonetheless, such a challenging observation must be 
pursued. Of the small number of brown dwarf candidates 
identified by precision Dopplcr surveys to date (around 
11), many have been shown to actually be stars (Halb- 
wachs et al. 2000), making bona fide brown dwarf com- 
panions rare, and well worth unambiguously identifying. 

6. CONCLUSIONS 

We present first results from the Anglo- Australian 
Planet Search, which is now the highest precision Doppler 
planet search in the southern hemisphere. These include 
a new member of the "51 Peg" -like class of planets around 
the F8 dwarf HD 179949, and a new brown dwarf candidate 
companion to the GO dwarf HD 164427. This programme 
has now been running on the AAT for over 2 years and 
has demonstrated velocity precisions reaching 3 ms -1 . We 
confidently expect many more discoveries in the years to 
come - particularly as it extends its sensitivity toward the 
Jupiter-like planet regime. 

The Anglo-Australian Planet Search team would like 
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of the AAO, Dr Brian Boyle, and the superb tech- 
nical support which has been received throughout the 
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R.Patterson, D.Stafford, F. Freeman, S.Lee, J.Pogson and 
G.Schafcr. We further acknowledge support by; the 
partners of the Anglo-Australian Telescope Agreement 
(CGT,HRAJ,AJP); NASA grant NAG5-8299 & NSF 
grant AST95-20443 (GWM); NASA grant NAG5-4445 
& NSF grant AST-9619418 (SSV); NSF grant AST- 
9988087 (RPB); NASA grants NCC5-96 & NCC5-228, 
NSF grant HRD-9706268 & the Richard Lounsbery Foun- 
dation (GWH); and Sun Microsystems. 
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Table 1 
Velocities for HD 179949 

JD RV Uncert 

(-2440000) (ms- 1 ) (ms" 1 ) 



11120.910 


118.1 


5.8 


11383.024 


-43.3 


8.9 


11410.943 


-47.7 


5.3 


11413.052 


-15.0 


4.7 


11413.924 


-42.1 


6.1 


11472.917 


-3.6 


5.7 


11683.133 


-40.6 


8.0 


11684.180 


113.7 


6.1 


11706.151 


91.1 


9.8 


11718.161 


127.1 


6.8 


11742.960 


115.9 


5.9 


11743.992 


-42.9 


5.9 


11745.101 


0.0 


5.8 


11766.970 


40.8 


6.2 


11768.000 


109.6 


8.5 


11770.107 


44.6 


6.1 


11792. 751 a 


122.2 


2.6 


11793.798 a 


-80.7 


2.5 


11795.803 a 


115.6 


2.5 


11827.956 


-74.3 


6.1 


11828.911 


61.5 


6.7 


11829.915 


117.1 


8.2 


11855.942 


-74.4 


10.3 



Observations from Keck. All 
other observations are from the 
A AT. 
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Table 2 



Orbital Parameters 



Parameter 


HD 179949 


HD 179949 


HD 164427 




(e=0.0) 


(e fitted) 




Orbital period P (d) 


3.093±0.001 


3.093±0.001 


108.55±0.04 


Velocity amp. K (ms _1 ) 


101.3±3.0 


102.2±3.0 


2229±77.0 


Eccentricity e 


0.0 (fixed) 


0.05±0.03 


0.55±0.02 


u>n 


0.0 


226±25 


356.9±0.5 


a\ sin i (km) 


4307.3±15.0 


4338.7±15.0 


(2.776±0.049)xl0 6 


Periastron Time (JD-244000) 


11718.19±0.1 


11723.24±0.1 


11724.6±0.2 


Msini (Mj UP ) 


0.84±0.05 


0.84±0.05 


46.4±3.4 


a (AU) 


0.045±0.004 


0.045±0.004 


0.46±0.05 


RMS about fit (ms" 1 ) 


10.8 


10.3 


7.7 
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Table 3 
Velocities for HD 164427 



JD RV Uncert 

(-2440000) (ms- 1 ) (ms" 1 ) 



10917 


.287 


-1071 


.2 


7.5 


11002 


.091 


-897. 


.7 


10.2 


11118 


.891 


-1026 


.5 


9.8 


11119 


.905 


-1025. 


.7 


8.3 


11236 


.281 


-1065 


.7 


6.4 


11274 


.300 


-3. 


9 


5.3 


11385 


.868 


347. 


,2 


6.8 


11386 


.879 


513. 


.8 


4.9 


11410. 


.911 


703. 


.9 


4.2 


11413 


.034 


378. 


,0 


3.6 


11413 


.904 


286. 


,7 


5.6 


11472 


.893 


-914. 


.8 


4.6 


11630 


.300 


391 


.3 


4.9 


11683 


.096 


-1014 


.7 


5.6 


11684 


.147 


-1006 


.3 


6.4 


11706 


.112 


-217. 


.9 


4.9 


11717. 


.923 


1718. 


.7 


5.3 


11718. 


.123 


1775. 


.8 


4.9 


11742 


.913 


-44. 


.6 


5.5 


11743 


.933 


-126. 


.6 


7.0 


11745 


.067 


-218. 


.7 


4.6 


11766 


.942 


-968. 


.3 


5.2 


11767. 


.992 


-985. 


.4 


8.4 


11827. 


.914 


2169. 


.0 


5.5 


11828 


.890 


2480. 


.8 


6.3 


11829. 


.894 


2791. 


,5 


7.0 


11855 


.921 


-371. 


.0 


8.0 



Tinney et al. 



Table 4 

Predicted K magnitudes for HD 164427 companion 



Companion Mass 


Age 


M K 


AM K 


HD 164427 M K =5A 






0.18M Q a 




7.7 


2.3 


0.10M Q a 




9.2 


3.8 


75Mjup b 


lGyr 


11.2 


5.9 




5Gyr 


11.7 


6.3 


42 Mjup b 


lGyr 


13.7 


8.3 




5Gyr 


>16 


>10.6 



a Empirical Mk values due to Henry & 
McCarthy (1993). 

b Model Mk values due to Chabrier et al. 
(2000). 



